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1 . INTRODUCTION 

The strong interaction between structural dynamics and active control is a 
well-recognized challenge in the analysis of controlled flexible structures. 
But it is only recently that the same interaction has been exploited in the 
design process. The traditional design approach in which the control design 
comes very late in the development — typically after the structure has been 
designed and built — is no longer viable. Although this approach has produced 
satisfactory results for the attitude control of relatively rigid space 
structures, it will not be suitable for the ambitious space missions that 
require precise controlled pointing of telescopes, interferometers and the 
vibration suppression for science instruments mounted on large flexible 
structures. In such systems, designing the structure and designing its control 
become entwined. This dictates early consideration of the control design — 
well before any detailed structural design is finalized. And just as the 
structure is optimally designed to meet such performance metrics as minimum 
mass or response to external disturbances, it should be optimally designed to 
meet its ultimate control performance as well. 

A natural way to introduce the control element into the overall design 
process is through an optimization procedure that combines the structural and 
control design criteria into a single problem formulation. A number of authors 
[1-6] have taken this perspective. In terms of the types of design parameters 
and constraints considered, Ref. (2) is probably the most extensive in that the 
design variables include structure parameters, actuator locations and feedback 
matrix. Static output feedback is used, and the performance objectives include 
total mass and robustness measures. Constraints are imposed on the eigenvalue 
placements, performance bounds, and structural constraints. Since not all of 
the constraints are commensurate, they are relaxed using a homotopy approach. 
Just as with Ref. (6), the approach taken in the present paper is not to 
produce the "best" optimized point design, but to produce a family of Pareto 
optimal designs representing options that assist in early trade studies. The 
philosophy is that these are candidate designs to be passed on for further 
consideration, and their function is more to guide the system design rather 
than to represent the ultimate design. 

An optimization approach that is consistent with this philosophy is to 
utilize multiple cost objectives that include an 1QG cost criterion in 
conjunction with structural cost(s), and possibly other control related costs. 
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After introducing the combined objective formulation in the setting of vector 
optimization, we derive the necessary conditions for Pareto optimality. No 

behavioral constraints are explicitly imposed in the problem formulation and a 
homotopy approach is used to generate a family of optimum designs. The intent 
of this paper is to explore this design approach and to provide an exposition 
of the computational aspects of the problem using numerical examples. 

2. COMBINED OPTIMIZATION 

The combined optimization approach can be best appreciated when contrasted 
with the traditional sequential optimization. In the sequential optimization, 
the structure is first optimized by selecting the design variable, a (e.g. 
member sizes) which minimize a structural criterion J s (a) - often taken as the 
mass of the structure subject to some behavioral constraints h(a)>0 on 
deformations, stresses, open-loop frequencies, etc. 

min J (a) ; h(a)>0, aeD (2.1) 

a 

where D is the physical domain for a. Second, having completely specified the 
optimal structural design a*, the control optimization is carried out with a* 
fixed. For example, LQG or optimal control designs pose the problem: 

min J (a*,C) (2.2) 

c c 

where J c represents either of the control criterion, and C is allowed to vary 
over the class of stabilizing compensators for the plant. 

By contrast, in the combined optimization formulation, the goal is to first 
merge the criteria of interest (here J s and J c ) into one using non-negative 
multipliers p, and 6, then optimize the combined criterion over the original 
design variable space a, C: 

min [0J (a) + 6J (a,C] (2.3) 

Q,C S C 

The following expression compares the results of the two optimization 
procedures outlined above. 

min [0J (a) + 6J (a,C)]<Fmin pj (a) + min 6J (a*,C)] (2.4) 

a , C S C a S C 

The right-hand side of (2.4) corresponds to performing the sequential 
optimization by solving (2.1) for a*, followed by solving (2.2) for C*. Note 
that the optimal solution of the right-hand side is independent of p and 6- 
but not so for the combined optimization embodied by the left-hand side. In 
terms of the total cost, expression (2.4) states the fact that the combined 
optimization is never inferior to the sequential optimization. In the vector 
optimization setting, the relative weighting of p and 6 serves as a parameter 
that allows the generation of an entire family of Pareto optima. 

In the present paper, only two objective criteria are dealt with. But it 
is not difficult to generalize the approach to incorporate other criteria such 
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as minimum open- loop frequency and certain controller robustness measures. In 
general these criteria are noncommensurate , and there is no unique solution 
that minimizes all criteria simultaneously. Thus, one must look for 
Pareto optimal solutions as outlined below. 


First one assembles the criteria J^:D->R, i~l,...,N into a single criterion 
J:EHR N , J(a) = (J]^(a) , . . . , J^(a)) T - Tlien the cone C Q = {xcR N :x^>0, i«l,...,N} is 
defined to induce a partial ordering < on R N by x<y if y-x eC G . Now let aeD. 
A design vector a*£D is said to be (strongly) Pareto optimal if J(a) < J(a*) 
implies J(a)=J(a*). A necessary condition for Pareto optimality is contained 
in the following theorem due to Lin [7]. 


Theorem 2.1: If a* is Pareto optimal for the combined criterion J, and D is an 
open set, then there exists a nonzero vector ZeC Q such that Z^ J * ( a* ) -0 . Here 
J* denotes the differential of J . 


For the two-term optimization problem in (2.3), we find that the Pareto 
optimal solution to J = (J s , J c )^ can be generated by solving for the necessary 
conditions for extremizing the following convex combination J^: 

J A - (1-A) J g (a) + AJ c (a,C); Ae[0,l] (2.5) 


where A replaces y3, 8 in (2.3). The form of Eq. (2.5) suggests a homotopy (or 
continuation) approach for generating a family of Pareto optima as A is 
propagated from 0 to 1 . 

3. NECESSARY CONDITIONS 


We begin with the n s degree -of- freedom dynamical system 

M(a)r + D(a)r + K(a)r - G^u + G 2 V (3.1) 


where M, D, and K are the n s xn s mass, damping and stiffness, G ^ is the constant 
n s xn u control influence matrix, and G 2 is the constant n s xn<i disturbance 
matrix. The vectors r, u, and v are respectively, physical degrees-of-freedora, 
control forces and disturbances. Let x - (r,r) . Then the first order state 
equation is 


x — A(a)x + B^(a)u + B2 (ck) v + u' 
where 



- 0 

I ' 


■ 0 

r 0 1 

A - 

-M" X K 

-M' X D 

. Bj. - 

M -1 Gi 

’ ® 2 ” M* 1 G 2 J 


(3.2) 


(3.2a) 


and an additional disturbance v' independent of a has been introduced for 
greater flexibility of the formulation. We assume that (3.2) has measured 
output variables y and controlled output variables z: 


y - C^x + w, 
z - C 2 X 


157 


(3.3) 



and that u, v * , and w are uncorrelated white noise processes with intensities 
Q v , V, and Q w , respectively. 

In the remainder of this paper, the total mass of the structure is assumed 
to represent the structural criterion J s . Thus, for a structure consisting of 
n a one -dimensional finite elements, each having a cross-sectional area a^, 
length and density p: 


J 

s 



p2 .a . 
r l l 


(3.4) 


For the control criterion J c associated with (3.2) and (3.3), we assume the LQG 
index 


J - lim E{z^(t)Dz(t) + x^D-x + u^(t)Ru(t)} 
c l 

t-*» 


(3.5) 


where E is the expectation, D and are non-negative definite weighting 

matrices, and R is positive definite. Although is assumed independent of a, 
D could possibly depend on a. The latter case would arise, for instance, if 
the first term in (3.5) were to represent the total energy in the system with 
z - (r,r)^ and D - diag(K,M). Under standard assumptions of stabilizability 
and detectability, the optimal compensator C* for (3.5) is implemented by [8]: 


u 

o 


T 

-B^Px 


o 


x - (A-K^C)x + Kj-y + B-u , 
o f o f y 1 o 


(3.6) 


and the optimal cost J c associated with this compensator is 


J (C) - tr{ P(B 0 Q + V) + P.PB 1 R' 1 B:F P). (3.7) 

c 2 v 2 ill 

where P and Pf are the unique positive definite solutions to the algebraic 
Riccati equations 


A T P + PA + D - PBtR^B^ P - 0, 
oil 

(3.8) 

AP f ♦ + B 2 Q^ ♦ V - p f c* q; 1 c 1 P f - 0 

and D q - cj DC 2 + and K f - P f C^ C^ 1 

With the above representation for J c , we seek the optimality conditions for 


min J (a)— (1- A) J (a)+<v ,a>+Atr{P(B 0 Q B^+V)+P-PB 1 R' 1 bJ : P) (3.9) 

A s z v z r i i 

a 
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subject to the constraints in (3.8). In (3.9) an n a - dimensional vector vc R 
has been introduced to "regularize" the problem and to serve the purpose of 
initializing the homotopy path. 

Proposition 3.1 Let E+ denote the set of n x xn x positive definite matrices. 
The optimality conditions for a* to be a local extremum of (3.9) subject 
to (3.8) are the zeros of the function H: RxR na x Z+ x E+ -+ R na x E+ x E+ 

defined by H - [H]^ } H 2 , H 3 , H 4 , H 5 ] , where 


H 1 (A,a f Z 1 ,Z 2 ,P,P f ) - 

3J 8(B,R'V) ,,T 3D 

d-Dv - 5 +v.+xtr( p — r*- 2 —* pp f p — \ — + z,t^-p + p“- + r - 2 
5a . 1 da, f da . 1 5a. da. da, 

1 i i i 1 i 


3<B R -1 B P ) 3A 3 * T 8 < B 20 „ B 2 > 8(C 1 CV 

• P t 1 P ' + Vt P f + P f 3.7 + 3l ' P f 3a. 1 P £l>. 

1 ill i 


i-1 n. 


(3.10a) 


H 2 (A, a, Z 1 ,Z 2 , P, P f ) - 

(A-BjR* 1 B^P) Z 1 +Z 1 ( A T - PBjR' 1 B^)+B 2 Q^B 2 +P f PBjR* 1 B^+B 1 R* 1 B^PP f (3.10b) 


, a, Z^ , Z 2 i P, P^-) 


(A 1 -^^ 1 C 1 P f ) Z 2 +Z 2 (A-P f C^Q w 1 C x )+ PB^^B^P 


(3.10c) 


H 4 (A, a, Z r Z 2 , P, P f ) 


a t p+pa+d -pb^^bTp, 

o 1 1 


H 5 (A, a. Z r Z 2 , P, P f ) 


AP f +P f Al+B 2 Q u B 2 +V - P f C I Q ; 1 C l P f 


(3 . lOd) 


(3 . lOe) 


The proof of proposition (3.1) is omitted here to conserve space, but is given 
in detail in Ref. [9]. Thus, for the LQG formulation, the necessary conditions 
involve solving two algebraic Riccati equations (3.10d, e) , two Lyapunov 
equations (3.10b, c) , and a gradient equation (3.10a) for a^, i— n a . In 
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the case of the LQR formulation, it is easily verified that the optimization 
statement expressed by (3.9) reduces to 


min J A (a)-(l-A)J s (a)-Ki/,a>+Atr{P(B 2 Q u B 2 +V) } (3.11) 
a 

and that the optimality conditions reduce to finding the zeros of the simpler 
function H-[H]^, H 2 , H 3 ] involving one Riccati and one Lyapunov equation 
(instead of two as in LQG) , in addition to the gradient equations: 


dj 3 < b 2 W 

H(A,a,Z,P) - (1-A) + u ± + Atr { P ^ ~ 


dD fl(B 1 R* 1 Bb 

+Z I 2 P fc: + i^- P 8 a t P »' 


i -1 n 


(3.12a) 


H 2 (A,a,Z,P) - A c Z + ZA T c + B^B* + V, 


(3.12b) 


H 3 (A,a,Z,P) - A T P + PA + D q - PB^’^P, 


(3.12c) 


with A c - A - BjR’^P. 


4 . HOMOTOPY STRATEGY 

For all Ae[0,l], our goal is to minimize (3.9) in the case of the LQG 
formulation or (3.11) in the case of the LQR formulation by finding the design 
variables a that satisfy the corresponding optimality conditions ( 3 . 10 ) or 
(3.12). The basic strategy is: given the solution at a value A 0 , smoothly 
propagate it to a new solution at A 0 +AA via some local mechanism such as Newton 
method or iterative optimization. This strategy has been analyzed in detail in 
Ref. (9). In the following, only a summary of the results is given without 
proofs, assuming the LQR formulation. 

Let x denote a generic point (a, Z, P)c R na x ][+ x £ + so that H(A, x) stands 
for H(A, a, Z, P) . In determining the zeros of H, the following proposition 
asserts that in a small neighborhood about the optimal at A-0, there is a 
smooth path parameterized by A consisting of the global optimal solution. 

Proposition 4.1: Suppose that min J s (a) has a unique global solution a* 
satisfying the second order sufficiency condition on the positivity of the 
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Hessian [J s (a*)], a a >0. Further, assume that J s is coercive so that 

i j 

|J s (a)|-^° as Then there exists €>0 such that (3.11) has a unique global 

solution for A<e . 

The next proposition provides a sufficient condition for the path to remain 
locally optimal. 

Proposition 4.2: Let A) - (A, x(A)) denote a smooth path in 

[0 , r)xR na x X+ x Z+ with H(<£(A)) - 0 and H, x invertible for Ae[0,r). Such an r 
is guaranteed by Proposition 4.1. Then a(A) is a local minimum for for each 
Ac [0 , r) . 

The purpose of the following lemma is to demonstrate that the zero set of H 
is "genetically” well-behaved. 

Lemma 4.3: Suppose that H(0,x)=0 has a unique solution. Then for almost every 

choice of (i/ , V, D^)cR na x £+ x £+, the solution path emanating from (0 x*) is 

diffeomorphic to the real R and every other component of H'-^o) is 
diffeomorphic to either a circle or R. 

Thus, following the path defined by one of the zero curves of H, not just 
the one emanating from the optimal at A-0, will not lead to a pathological 
behavior such as bifurcations or curves with infinite length in bounded sets. 
Another fundamental and generally difficult question that arises when employing 
homotopy methods is whether or not the path remains bounded. The following 
result provides a partial answer to this question. 

Proposition 4.4: Suppose that J s , B^, 0 o , and A are all polynomials in 

ai,...a na » and assume coercivity of J s (a)/|a|. If H(0,x)-0 has a unique 
solution, then for any e>0 and for almost every triple (i/, V, D^) 
cR na x x £+, the component of H“^(0) containing (0,x*) is a bounded set in 
[ 0 , 1 - c ] xR na x{ + x X+. 

5. NUMERICAL RESULTS 

The numerical experiments described in this section demonstrate the results 
of the foregoing theory. Two prototype examples are used; both employing the 
LQR formulation. Implementation of the homotopy strategy of the previous 
section is achieved by iterative optimization, wherein the solution path for 
minimizing (3.11) in terms of the homotopy parameter A starts at A— 0 with 
,a na initialized to a predetermined sufficiently small allowable size a c . 
At this point in the solution space, is fully weighted toward minimizing J s 
only. Minimization of J^ Q thus yields a G * for which H(A 0 , x o *)-0. For the 
next iteration and for every succeeding one, A is incremented and the weighting 
is shifted gradually toward J c . For a typical iteration j, the following steps 
are performed: 

(i) Aj * A j -1 + AA j -1 

(ii) Initialize the minimization of J^j by using This will 

result in aj* for which conditions (3.12) hold. 
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In performing the minimization in (ii) above, we employed the Automated Design 
Synthesis (ADS) system of general purpose subroutines [10] . ADS provides a 
wide selection of options at three levels: strategy, optimization, and line 

search. Available strategies include sequential linear and quadratic 
programming, and sequential unconstrained minimization coupled with various 
penalty methods. At the optimization level, one can choose between Fletcher - 
Reeves algorithm and the variable metric method of Broydon- Fletcher -Goldfarb- 
Shanno (BFGS) for unconstrained minimization, or Zoutendijk's method of 
feasible directions and modifications thereof for constrained minimization. 
For one -dimensional line search, the options include a combination of 
polynomial interpolation/extrapolation, solution bounding, and the method of 
Golden Section. Not all combination of options are compatible at the three 
levels, and the program parameters must be adjusted to suit the problem at 
hand. For this purpose, an analytical function was contrived which possessed 
such features as: easy to compute closed form solution, multiple minima, and 

insensitivity of the functions gradient near the minima to design parameters. 
Several compatible options were tried and the program parameter values (e.g., 
move limits and convergence criteria on the absolute and relative changes in 
objective function between two successive iterations) were adjusted until the 
closed form and numerical solution agreed within as few iterations as possible. 
As a result of these numerical experiments, the popular BFGS variable metric 
method for unconstrained minimization emerged as the one of choice for use in 
the combined control -structure optimization examples that follow. During the 
one dimensional line search, the minimum is located by first computing the 
bounds, then using polynomial interpolation. 

Example 1: The cantilever beam of Fig. 1 resembling a flexible appendage* of a 

large structure is modelled by three finite elements with two degrees- of - 
freedom (dof) at each node. The structural design variables are the x- 
sectional areas , 02 , 03 of the elements. The disturbance v represents a 
transverse sound pressure modelled by uncorrelated unit impulses at t -0 
concentrated at the three nodal transverse dof. The control force u is applied 
at the free end along the transverse dof direction. With the J s given by 
(3.4), we seek the minimum of (3.11) for Ae[0, 1] , subject to conditions (3.12). 
Here, the weighting matrices were taken as D- 10 ^ x diag(K,M), and R-10’\ and 
the initial a ^ used were a^-a o -lxl0“^, i~l,...,3. 

Table 1 lists the family of Pareto optimal designs cr^* that minimize 
J^, Ae[0,.99] along with the corresponding values for J c , J s and J^. The 

variations of the Pareto optimal J C *(A), J S *(A) and J^*(A) are shown in Fig. 2. 
A number of observations can be made from Table 1 and Fig. 2: 

(i) The noncommensurate nature of the two costs J c and J s is apparent as the 

weight is shifted between them: while J c is a strictly decreasing 

function of A, J s is a strictly increasing function of A. This is 
consistent since a stiffer structure requires less control energy. 

(ii) Except near A-0, the optimal structural shapes that minimize for the 

disturbance, choice of D and R, and control forces described above are 
essentially near the fixed end, and a much smaller 03 near the 

free end. This is a physically plausible optimum shape that minimizes 
the mixture of high strain energy density near the fixed end and high 
kinetic energy density near the free end. Other choices of disturbance, 
control location and D, R are expected to alter the optimal beam shape. 
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(iii) Although the design at A*0 is guaranteed to be globally optimal 
(Proposition 4.1), it is possible that designs generated as A is 
continued may be only locally optimal. To assess this possibility, the 
optimal designs listed for A-.200, A-.700 and A-.900 were re-examined 

separately. For each case, the minimization was restarted with randomly 
selected initial values. In most cases, the minimization converged 
to the same or to a higher minimum than obtained in Table. 1. 


Table 1. Pareto Optimal Designs of Example 1 


A 

Optimal Design xlO"^ 

«1 Q 2 “3 

^c 

Js 


0.0 

.00010 

.00010 

.00010 

350000. 

.007 

.007 

. 000001 

.00035 

.00038 

.00032 

96000. 

.026 

.122 

.00001 

.00116 

.00114 

.00077 

30800. 

.076 

.384 

.0001 

.00570 

.00625 

.00287 

5870. 

.369 

.955 

.005 

.0585 

.0708 

.0187 

495. 

3.68 

6.14 

.010 

.0801 

.0802 

.0212 

391. 

4.52 

8.89 

.020 

.113 

.111 

.0300 

265. 

6.33 

11.52 

.040 

.159 

.155 

.0438 

177. 

8.90 

15.64 

.100 

.237 

.249 

.0761 

102. 

14.00 

22.77 

.200 

.367 

.349 

.112 

61.8 

20.62 

28.86 

.300 

.425 

.457 

.154 

46.1 

25.80 

31.92 

.400 

.563 

.532 

.189 

35.50 

31.97 

32.99 

.500 

.622 

.650 

.231 

27.93 

37.43 

32.68 

.600 

.795 

.751 

.279 

21.33 

45.45 

30.97 

.700 

.886 

.934 

.346 

16.86 

53.92 

27.98 

.800 

1.20 

1.14 

.44 

11.78 

69.39 

23.30 

.900 

1.68 

1.59 

.62 

7.30 

96.98 

16.27 

.940 

1.93 

2.07 

.82 

5.40 

119.9 

12.27 

.980 

3.37 

3.22 

1.37 

2.58 

198.2 

6.49 

.990 

3.95 

4.40 

1.81 

1.81 

253.6 

4.33 


Example 2: The beam in this example (See Fig. 3) simulates a flexible 
appendage (length * 45 m) attached to a rigid hub (inertia - 50 kg.m 2 ) to which 
a control torque is applied to counteract the transverse unit impulse at the 
free end. The beam is modelled by three finite elements of constant width — 
.001m, but whose nodal depths dj L ,...,d 4 represent design variables having a 
lower bound ** .001 m. Here again, J s represents the total mass of the flexible 
beam (excluding the hub). For the control objective J c , the response energy is 
weighed by so as to minimize the transverse free end displacement, and R is 
taken - l.xl0“V 

Table 2 and Fig. 4 represent analogous results to those presented for 
Example 1 in Table 1 and Fig. 2. In addition to observation (i) of the 
previous example - which holds here as well - the following remarks can be made 
with reference to the results of this example: 
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(i) For small values of A (e.g. A~0.1), where the total mass J s dominates 

the minimization of J^, the optimal shapes tend to have a small slope 
from d]^d 2 “ > d 3 , followed by a sharper slope from d 3 ~^d 4 . As A increases, 
minimizing J s becomes less important than minimizing J c (tip 
displacement response energy plus control energy). As a result, the 
beam becomes gradually stiffer, and the monotonic slope from 
associated with small A values gradually disappears at A=.45, then gives 
way to a pronounced inflection of slope for for A>.45. This 

results in a larger allocation of mass at the tip. This type of shape 

is physically consistent with the requirements of the two parts of the 
control objective J c : a stiffer structure near the hub that is reduced 

toward the tip (free end) makes best distribution of mass, while 
minimizing the tip displacement response. On the other hand, a large 
mass at the tip (where the disturbance exists) makes the disturbance 
less effective - thus requiring less control effort. 

(ii) To confirm the above interpretation, the case of A=0.7 in Table 2 (i.e. 

R-10"^) was re-examined for smaller and larger values of R; R-10"° and 

R-10"^, respectively. As Fig. 5 shows, smaller values of R give more 
weighting to the tip displacement response energy part of J c , thus 
giving rise to the optimum shape being a stiffer structure near the hub 
which is reduced toward the tip. Conversely, larger values of R (e.g. 
R-10~2) give more relative weighting to the control energy part of J c , 
which is best minimized by the presence of the heavier tip mass. It is 
interesting to note the similar effect of varying R and varying A on the 
optimal shapes. 


Table 2. Pareto Optimal Designs for Example 2. 


A 

d l 

Optimum Design 
d2 d 3 

d 4 

Jc 

Js 

J A 

.000 

.001 

.001 

.001 

.001 

1 . 6xl0 +10 

.075 

.075 

.0001 

.02404 

.02363 

.02291 

.01366 

3355.26 

1.628 

1.963 

.001 

.03552 

.03485 

.03359 

.02065 

482.6 

2.404 

2.884 


.05223 

.05134 

.04940 

.03043 

70.78 

3.54 

4.21 

.100 

.07699 

.07559 

.07157 

.05303 

10.03 

5.28 

5.757 

.200 

.08759 

.08563 

.08031 

.06659 

5.33 

6.05 

5.906 

.300 

.09550 

.09267 

.08657 

.07847 

3.54 

6.62 

5.703 

.400 

.10258 

.09867 

.09187 

.09085 

2.56 

7.15 

5.31 

.500 

. 10944 

.10389 

.09670 

.10470 

1.93 

7.66 

4.79 

.600 

.11715 

.10948 

.10135 

.12176 

1.46 

8.22 

4.17 

.700 

.12657 

.11576 

.10681 

.14516 

1.08 

8.92 

3.43 

.800 

.13944 

.12322 

.11304 

.18055 

.77 

9.87 

2.59 

.900 

.16369 

.13641 

.12430 

.24923 

.47 

11.63 

1.59 

.920 

.17299 

. 14104 

.12893 

.27374 

.41 

12.28 

1.355 

.940 

.18542 

.14815 

.13565 

. 30544 

.34 

13.18 

1.107 

.960 

.20578 

.16117 

.14818 

.35131 

.26 

14.64 

0.83 

.980 

.25035 

.19468 

.18421 

.42401 

.16 

17.83 

0.52 

.990 

.29401 

.25488 

.24878 

.48165 

.09 

22.20 

0.32 
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(iii) Of general interest to problems in combined optimization - at least 
numerically - is the question as to the degree of "roughness" of the 
hyper-surface J^(a) . A partial answer to this question is provided in 
Fig. 6 after introducing idealizations that reduce the number of 

variables from four (d^, . . . t d^) to only two (d 3 , d 4 ) , so that a three 
dimensional plot could be generated. Figure 6 shows such a surface in 
the neighborhood of the optimum for the case A— 0.7 in Table 2. This is 
achieved by fixing d^-.13, allowing d 3 and d 4 to assume various values 
larger and smaller than those in Table 2 for A-0.7, and letting 
d 2 “h(d]^+d 3 ) . Assuming that the idealizations above (which led to 

reducing the dimensionality of the surface) did not alter the basic 

topology of surface, it appears from Fig. 6 that is a smooth 

function of the design variables - at least in the neighborhood of the 
minimum. Furthermore, with these idealizations it appears that is 

relatively flat near the minimum along the d 4 axis, and that the optimum 
shape is some linear combination of the four basic shapes depicted at 
the corners . 

6 . CONCLUSIONS 

An approach for combined control -structure optimization keyed to enhancing 
early design trade-offs has been outlined and illustrated by numerical 
examples. The approach employs a homotopic strategy and appears to be 
effective for generating families of designs that can be used in these early 
trade studies. 

Analytical results were obtained for classes of structure/control 
objectives with LQG and LQR costs. For these, we have demonstrated that global 
optima can be computed for small values of the homotopy parameter. Conditions 
for local optima along the homotopy path were also given. Details of two 
numerical examples employing the LQR control cost were given showing variations 
of the optimal design variables along the homotopy path. The results of the 
second example suggest that introducing a second homotopy parameter relating 
the two parts of the control index in the LQG/LQR formulation might serve to 
enlarge the family of Pareto optima, but its effect on modifying the optimal 
structural shapes may be analogous to the original parameter A. 
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Figure 1 


EXAMPLE 1. 

CANTILEVER BEAM PROBLEM 



STRUCTURAL MODEL: MASS DENSITY « 1660 Kg/m*. MODULUS * 9.56 x 10 10 Him 1 
MODAL DAMPING * 0.5% 

CONTROL: DISTURBANCE a TRANSVERSE PRESSURE IMPULSE CONCENTRATED 

AT THE NODES 

RESPONSE ENERGY WEIGHTED BY 5 > Olag (K.M) x 10 2 , 
CONTROL ENERGY WEIGHTED BY R « 1 x 10"* 

DESIGN VARIABLES: a,, a„ a, * 1 x 10 7 

Figure 2 

CANTILEVER BEAM OPTIMIZATION 
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Figure 3 


EXAMPLE 2. 
HUB-BEAM PROBLEM 



STRUCTURAL MODEL: MASS DENSITY . 1660 Kg/m 3 , MODULUS » 9.56 x 10 10 N/m 2 
MODAL DAMPING > 04%, 

CONTROL: DISTURBANCE * UNIT IMPULSE 

RESPONSE ENERGY WEIGHTED TO MINIMIZE END DISPLACEMENT, R > 1 x 10' 4 

DESIGN VARIABLES: d v * • •, d 4 20.001 


Figure 4 

HUB-BEAM OPTIMIZATION 
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